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We discuss the behavior of quantum and classical pairwise correlations in critical systems, with 
the quantumness of the correlations measured by the quantum discord. We analytically derive 
these correlations for general real density matrices displaying Z2 symmetry. As an illustration, we 
analyze both the XXZ and the transverse field Ising models. Finite-size as well as infinite chains are 
investigated and the quantum criticality is discussed. Moreover, we identify the spin functions that 
govern the correlations. As a further example, we also consider correlations in the Hartree-Fock 
ground state of the Lipkin-Meshkov-Glick model. It is then shown that both classical correlation 
and quantum discord exhibit signatures of the quantum phase transitions. 
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I. INTRODUCTION 

The concept of correlation, i.e., information of one 
system about another, is a key element in many-body 
physics. Indeed, the properties of a many-body sys- 
tem are strongly affected by changes in the correlations 
among its constituents. These changes are responsible for 
the occurrence of remarkable phenomena such as a quan- 
tum phase transition (QPT), which is a critical change in 
the ground state of a quantum system due to level cross- 
ings in its energy spectrum. QPTs occur at low temper- 
atures T where the de Broglie wavelength is greater than 
the correlation length of the thermal fluctuations (effec- 
tively T = 0) [lj, I2I . Striking examples of QPTs include 
metal-insulator transitions in strongly correlated elec- 
tronic materials, magnetic transitions in quantum spin 
lattices, superfluid-Mott insulator transitions in atomic 
gases induced by a Bose-Einstein condensation, among 
others. 

Correlations can be both from classical and quantum 
sources. The existence of genuinely quantum correlations 
can be usually inferred by the presence of entanglement 
among parts of a system. Indeed, entanglement displays 
a rather interesting behavior at QPTs [3(, being able to 
indicate a quantum critical point (QCP) through nonan- 
alyticities inherited from the ground state energy [J, || . 
Moreover, for one-dimensional critical systems, ground 
state entanglement entropy exhibits a universal logarith- 
mic scaling governed by the central charge of the Vi- 
rasoro Algebra associated with the underlying confor- 
mal field theory [|| 0, H| . Remarkably, the logarithmic 
scaling is robust against disorder [{J, [Tfl El, EH- ^ n 
higher dimensions, entanglement usually scales following 
an area law for noncritical systems (see, e.g., Ref. [l3j). 
For critical models, violations of the area law have been 
found [3, EH , with logarithmic- type corrections appear- 
ing. More recently, it has been observed that an area law 
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is generally implied by a finite correlation length when 
measured in terms of the mutual information [16j |. This 
remarkable behavior of entanglement at criticality is in- 
deed a consequence of the correlation pattern exhibited 
by the ground state of the system. 

Nevertheless, although entanglement provides a route 
to find out the existence of quantum correlations, it can 
be shown that they can appear even when entanglement 
is absent [l7], EH- Quantum correlations, which can be 
measured by the quantum discord [l7j . often arise as a 
consequence of coherence in a quantum system, being 
present even for separable states. Moreover, entangled 
states commonly involve more than only quantum cor- 
relations, i.e., they usually carry classical correlations 
among their parts. In a multiparty mixed-state scenario, 
the possibility of nonseparable states with purely quan- 
tum correlations (with no supporting background of clas- 
sical correlations) have also been investigated but 
genuine indicators of multiparticle classical correlations 
are still under debate [2(|. By focusing on condensed- 
matter systems, the aim of this work is, starting from 
an arbitrary ^-symmetric model, to investigate pair- 
wise correlations, explicitly splitting up their classical 
and quantum contributions and analyzing their behav- 
ior at QPTs. As an illustration, ground states of both 
XXZ and transverse field Ising spin chains will be shown 
to share quantum as well as classical correlations for 
nearest-neighbor spin pairs, with both of them signalling 
the critical behavior of the system. An analysis of such 
correlations for the thermodynamic limit of XXZ and 
Ising chains has recently appeared in Ref. [2l[. Here, 
we generalize this analysis by providing analytical ex- 
pressions for the correlations in each model and also by 
discussing their finite-size behavior. As a further illus- 
tration, we describe the correlations in the Hartree-Fock 
ground state of many-body systems, taking the Lipkin- 
Meshkov-Glick (LMG) model as an example. 

The paper is organized as follows. In Sec.|H]we review 
the evaluation of correlations from the point of view of in- 
formation theory, describing in particular the concept of 
quantum discord. In Sec. lIIIl we discuss the computation 
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of classical and quantum correlations for general models 
displaying Z 2 symmetry. Sec. lIVi Sec. El and Sec.lVIIare 
devoted to illustrating our results for the XXZ chain, the 
transverse field Ising chain, and the LMG model, respec- 
tively. Conclusions are then presented in Sec. IVIII 

II. CLASSICAL CORRELATIONS AND 
QUANTUM DISCORD 

In classical information theory, the information ob- 
tained, on average, after knowing the value of a random 
variable X, which takes values within a set of proba- 
bilities {px}, can be quantified by its Shannon entropy 
H(X) = — J2 X logPz- We use the symbol log as denot- 
ing logarithm at base 2 throughout the text. By taking 
two such random variables X and Y , we can measure the 
correlation between them by their mutual information 

1(X:Y)=H(X) + H(Y)-H(X,Y), (1) 

where H(X, Y) = - J2x, y Pxy \°gPxy is the joint entropy 
for X and Y. By introducing the conditional entropy 

H(X\Y)=H(X,Y)-H(Y), (2) 

which quantifies the ignorance (on average) about the 
value of X given Y is known, we can rewrite Eq. (TT]) as 

X(X :Y) = H(X) - H(X\Y). (3) 

In order to generalize the above equations to the quantum 
domain, we replace classical probability distributions by 
density matrices. Denoting by p the density matrix of 
a composite system AB and by p A and p B the density 
matrices of parts A and B, respectively, the quantum 
mutual information can be defined as 

I(p A :p B ) = S(p A )-S(p A \p B ), (4) 

where S{p A ) — — r £ip A \ogp A is the von Neumann en- 
tropy for subsystem A and 

S{p A \ P B ) = S(p) - S(p B ) (5) 

is a quantum generalization of the conditional entropy for 
A and B. A remarkable observation realized in Ref. [l7| 
is that the conditional entropy can be introduced by a 
different approach which, although classically equivalent 
to Eq. @ , yields a result in the quantum case that differs 
from Eq. ([5]). Indeed, let us consider a measurement 
performed locally only on part B. This measurement 
can be described by a set of projectors {B k }. The state 
of the quantum system, conditioned on the measurement 
of the outcome labelled by k, becomes 

Pk = —{I®B k )p(I®B k ), (6) 
Pk 

where p k = Tr[(J £g> B k )p(I ® B k )] denotes the probabil- 
ity of obtaining the outcome k and I denotes the identity 



operator for the subsystem A. The conditional density 
operator given by Eq. © allows for the following alter- 
native definition of the quantum conditional entropy: 

S(p\{B k }) = J2PkS(p k ). (7) 

k 

Therefore, following Eq. ([3]), the quantum mutual infor- 
mation can also be alternatively defined by 

J(p:{B k }) = S(p A )-S(p\{B k }). (8) 

Eqs. (H]) and © are classically equivalent but they are 
different in the quantum case. The difference between 
them is due to quantum effects on the correlation be- 
tween parts A and B and provides a measure for the 
quantumness of the correlation, which has been called 
quantum discord [13] ■ In fact, following Refs. [13, H2], we 
can define the classical correlation between parts A and 
B as 

C(p)=maxJ(p:{B k }), (9) 

{B k } 

with the quantum correlation accounted by the quantum 
discord, which is then given by 

Q(p) = I(p A : p B ) - C(p). (10) 

III. PAIRWISE CORRELATIONS FOR 
^-SYMMETRIC QUANTUM SPIN LATTICES 

We will consider here an interacting pair of spins- 1/2 in 
a spin lattice, which is governed by a Hamiltonian H that 
is both real and exhibits Z2 symmetry, i.e. invariance un- 
der 7r-rotation around a given spin axis. By taking this 
spin axis as the z direction, this implies the commuta- 
tion of H with the parity operator &)f =1 o-f, where N 
denotes the total number of spins and af is the Pauli 
operator along the z-axis at site i. Note that a number 
of spin models are enclosed within these requirements as, 
for instance, the XXZ spin chain and the transverse field 
Ising model. Disregardi ng spon taneous symmetry break- 
ing (see, e.g., Refs. [H, I24II25I] for a treatment of spon- 
taneously broken ground states), the two-spin reduced 
density matrix at sites labelled by i and j in the basis 
{| TT); I TT), I TT), I TT)}, with | t) and | |) denoting the 
eigenstates of a 3 , will be given by 

/ a / \ 

61 z . . 

p= z b 2 • ( H) 
V/ dj 

In terms of spin correlation functions, these elements can 
be written as 

a = \ (1 + Gj + Gi + G%) , 
bi = \{l + G\ - Gi - G%) , 
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b 2 = 


\ (1 - G\ + G{ 


~ ^Zz) ! 




d = 


~ (1 - G\ - G{ 


+ G%) , 




z = 








f = 


4 \PxX ~ Gyy) , 




(12) 



where G k z = (a k ) (k — i,j) is the magnetization den- 
sity at site k and G 1 ^ = (o % a o 2 ^) (a, f3 — x,y, z) denote 
two-point spin-spin functions at sites i and j, with the 
expectation value taken over the quantum state of the 
system. Note that, in case of translation invariance, we 
will have that G k z = G k (VA;, k') and, therefore, &i = b 2 . 
Moreover, observe also that the density operator given in 
Eq. (fTTj) can be decomposed as 



I ® I ■ 
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\Ci<7 



C4I ( 



C 5 <7 



(13) 



with 



In particular, for translation invariant systems, we have 
that c 4 = C5. In order to determine classical and quan- 
tum correlations, we first evaluate the mutual informa- 
tion as given by Eq. ^ . The eigenvalues of p read 



A 


1 
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(1 


+ c 3 ) 
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(1 


+ c 3 ) 


A 2 
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4 


(1 


-c 3 ) 


A3 


1 
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(1 


-c 3 ) 



(c 4 


+ c 5 ) 2 + ( Cl 


- c 2 ) 


(C4 


+ c 5 ) 2 + ( Cl 


- c 2 ) 


(C4 


- c 5 ) 2 + (ci 


+ c 2 ) 2 


(c 4 


- c 5 ) 2 + (ci 


+ c 2 ) 2 



(15) 



Therefore, the mutual information is given by 

3 

I(p) = S(p A ) + S( P B ) + Y,^\og\ a , (16) 



Ct=0 



where 



S(P B ) 



(rf log rf + rf log rf 



rf logrf 



log r. 



(17) 



(l+c 5 )/2, rf = (l-c 5 )/2, rf - (l+c 4 )/2, and 
rj" = (1 — c 4 )/2. Classical correlations can be obtained by 
following a procedure that is similar to those of Refs. (l8l 



with rf 
1 2 



l2l| , but applying it now for the case of the general density 
matrix given by Eq. (fTTj) . We first introduce a set of 
projectors for a local measurement on part B given by 
{B k = VUkV" 1 }, where {II fc = \k){k\ : k = 0, 1} is the set 
of projectors on the computational basis (|0) = | j) and 
|1> = I J.)) and V £ U{2). Note that the projectors B k 
represent therefore an arbitrary local measurement on B. 
We parametrize V as 



V = 



COS ; 



(18) 



where < 8 < ir and < 4> < 2ir. Note that 8 and 
</> can be interpreted as the azimuthal and polar angles, 
respectively, of a qubit over the Bloch sphere. By using 
Eq. ([6|) and the equation IlfcerTIfc = l) fc IIfc , with ^3 
denoting the Kronecker symbol, we can show that the 
state of the system after measurement {B^} will change 
to one of the states 



1 

Po = 2 



Cl 


= 2z 


f 2/, 






Pi = 


c 2 


= 22 


-2/, 








C3 


= a + 


d-bi 


-b 2 , 




where 


Ci 


= a — 


d-bx 


+ b 2 , 






C5 


= a — 


d + b x 


-b 2 . 


(14) 


qui 



I + J2lo^ j J <g> (VTIoV*) , (19) 
/ + ^<Zi j( 7 J ] ^(VTIiyt), (20) 



<7fc 2 = 
<7fc3 = 

with k = 0, 1 and 



(-l) fe ci 
(-l) fe c 2 

(-1)* 



1 + (-l) fc C 4 W3 
W 2 



1 + (-l) fe c 4 w 3 
C3W3 + (~l) fc c 5 



1 



-l) fe c 4 ui 3 



(21) 



wi = sin w cos < 
w 2 = sin sin q 
W3 = cos 6*. 



(22) 



Then, by evaluating von Neumann entropy from 
Eqs. CEU) and flU} and using that S'(VTIfcFt) = 0, we 
obtain 

S( Pk ) = — log — — log — ^— , 

(23) 



with 



\ 3=1 



(24) 



Therefore, the classical correlation for the spin pair at 
sites i and j will be given by 



C(p) — max S(p 
{B k } 1 



A\ (So + Si) (So — Si) 
C4W3 



(25) 
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where Sk — S(pk)- For some cases, the maximization in 
Eq. (j2"5)) can be worked out and an expression purely in 
terms of the spin correlation functions can be obtained 
(e.g., the XXZ and Ising chains below). In general, how- 
ever, C(p) has to be numerically evaluated by optimizing 
over the angles 9 and <j>. Once classical correlation is ob- 
tained, insertion of Eqs. (flUl) and (f2l>)) into Eq. (jTUJ) can 
be used to determine the quantum discord. 



IV. THE XXZ SPIN CHAIN 



Al = - (1 - C\ + C2 + c 3 ) , 

A 2 = | (1 + c\ - c 2 + c 3 ) , 
A 3 = J (1 + ci + c 2 - c 3 ) . 



(30) 



In order to compute C(p) and Q(/?) we write ci, c 2 , and 
c 3 in terms of the ground state energy density. By us- 
ing the Hellmann-Feynman theorem [28|, l29( for the XXZ 
Hamiltonian (12"6T), we obtain 



Let us illustrate the discussion of classical and quan- 
tum correlations between two spins by considering the 
XXZ spin chain, whose Hamiltonian is given by 



H 



XXZ 



J 

"o22 (°f°f+i 



Aafaf 



(26) 



where periodic boundary conditions are assumed, ensur- 
ing therefore translation symmetry. We will set the en- 
ergy scale such that J = 1 and will be interested in a 
nearest- neighbor spin pair at sites i and Concerning 
its symmetries, the XXZ chain exhibits t/(l) invariance, 
namely, cr*] = 0, which provides a stronger con- 

straint over the elements of the density matrix than the 
Zi symmetry. Indeed, U(l) invariance ensures that the 
element / of the reduced density matrix given by Eq. (JTTJ) 
vanishes. Moreover, the ground state has magnetization 
density G z = (<rj) =0 (Vfc), which implies that 



h = 



d= \{l + G zz ) , 

b2 = \ (1 - G zz ) , 
1 



Z — (G XX + Gyy) , 

/ = 0. 



(27) 



where, due to translation invariance, we write G Q /3 — 
( a a a ' l f3~ 1 ) (Y*)- Due to the fact that a — d, we will have 
that C4 = C5 = 0, which considerably simplifies the com- 
putation of classical and quantum correlations. More- 
over, we will have that c\ = c 2 = 2z and c 3 = 4a — 1. 
Then, the maximization procedure in Eq. (|25j) can be 
analytically worked out [18| , yielding 

C(p) = ^^log(l-c) + il±^log(l + C ), (28) 



Cl — C 2 — - (G XX + Gyy) 



de 

-^r-. e x 



dA 



c 3 = G zz = -2 



9A 



(31) 



where is the ground state energy density 



(ipo\H X xz\i'o) 



1 



(G, 



G V 



AG 2Z ), 

(32) 



with IV'o) denoting the ground state of Hxxz- Eqs. (j3T]) 
and (|32p hold for a chain with an arbitrary number of 
sites, allowing the discussion of correlations either for 
finite or infinite chains. Indeed, ground state energy as 
well as its derivatives can be exactly determined by Bcthc 
Ansatz technique [26| , which allows us to obtain the cor- 
relation functions c±, c 2 , and c 3 . In Fig. [TJ we plot clas- 
sical and quantum correlations between nearest-neighbor 
pairs for an infinite XXZ spin chain. 




FIG. 1: (Color online) Quantum and classical correlations for 
nearest-neighbor spins in the XXZ chain for L — > oo. 



with c = max (|ci|, |c 2 |, |c 3 |). For the mutual information 
I(p) we obtain 



where 



I{ P ) = 2 + J2 A, log Ai, 



Ao = | (1 - ci - c 2 - c 3 ) , 



(29) 



Note that, in the classical Ising limit A — * oo, we have 
a fully polarized ferromagnet. The ground state is then 
a doublet given by the vectors | TT " " T) an d | 41 ■ ■ ■ 4), 
yielding the mixed state 



TXTT 



T 1 + ^144 



4X44---4I- (33) 



Indeed, this is simply a classical probability mixing, with 
C(p) = I(p) = 1 and Q{p) — 0. The same applies for the 
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antiferromagnetic Ising limit A — > — oo, where a doubly 
degenerate ground state arises. Moreover, observe that 
the classical (quantum) correlation is a minimum (max- 
imum) at the infinite order QCP A = — 1. On the other 
hand, both correlations are discontinuous at the first- 
order QCP A = 1. This is indeed in agreement with the 
usual behavior of entanglement both at infinite and first- 
order QPTs. For an infinite-order QCP, entanglement 
commonly display a maximum at the QCP 3^, 3ll. |32|. 
while for a first-ord er Q CP, entanglement usually exhibits 
a jump at the QCP [UEsH]. Nevertheless, we note that in 
the specific case of the ferromagnetic QCP A = f and for 
pairwise entan glem ent measures such as concurrence (35j 
and negativity [36j, no j ump is detected, being hidden by 
the operation max [37]]. It is interesting to observe the 
behavior of the functions |ci| = \c2\ and | C3 1 that govern 
the classical and quantum correlations. For A < —I, we 
have that \c\\ — |ca| < |ca] , which means that the classi- 
cal correlation is governed by | C3 1 . For — I < A < 1, we 
have that |ci| = |ca| > | C3 1 , with the crossing occurring 
exactly at the infinite-order QCP. Therefore, the corre- 
lations are governed by different parameters in different 
phases. For A > 1, we obtain |ci| = |ca| = and |c 3 | = 1, 
which implies that C(p) — 1 and Q(p) = 0. These results 
are shown in Fig. [2] below. 




FIG. 2: (Color online) Values of the parameters |ci| = [ca | 
and I C3 1 as a function of A. The parameter |c3 1 governs the 
correlations for A < — 1 while |ci| = |ca| is dominant in the 
gapless disordered phase — 1 < A < 1. In the ferromagnetic 
phase, |ci| = |ca| = and |ca| = 1, which implies the fixed 
value C(p) — 1 and Q(p) = for any A > 1. 



with periodic boundary conditions assumed, namely, 
<t£ +1 = <rf. As before, we will set the energy scale such 
that J = I and will be interested in a nearest-neighbor 
spin pair at sites i and i + 1. This Hamiltonian is Z2- 
symmetric and can be exactly diagonalized by mapping 
it to a spinless free fermion model with single orbitals. 
This is implemented through the Jordan- Wigner trans- 
formation 

of = 1-2$*, 

of = -nl 1 - 2 ^) ( c *+ c i), ( 35 ) 

j<i 

where cj and Ci are the creation and annihilation fermion 
operators at site i, respectively. By rewriting Eq. ([34f in 
terms of cj and a we obtain 



1^ 

Hi = -J^2(4<k+i + 4+1* + 44+1 + 



-■"£( 



1 - 2c] c, 



(36) 



In order to diagonalize Hj we consider fermions in mo- 
mentum space 



c,e 



I L 



4 = 71^/'' 

.7=1 



-ikr-j 



(37) 



where ct and c k are creation and annihilation fermion 
operators with momentum k, respectively, and Vj is the 
fermion position at site j. The wave vectors k satisfy 
the relation ka = 2irq/L, where a denotes the distance 
between two nearest-neighbor sites and q = —M, —M + 
M-l, M, with M = (L-l) /2 and L taken, for sim- 
plicity, as an even number. Then, by inverting Eq. (|37[) 
and inserting the result in Eq. (|36p . we obtain 



(38) 



Hj = J ^ |^2 (ff — cos ka) CfcCfe 



-i sin ka 



+ C-kCk 



V. THE TRANSVERSE FIELD ISING MODEL 

Let us consider now the Ising chain in a transverse 
magnetic field, whose Hamiltonian is given by 



Diagonalization is then obtained by eliminating the terms 
c lfe c fe an d c_^Cfe from the Hamiltonian given by Eq. (|38[) . 
which do not conserve the particle number. This is in- 
deed achieved through the Bogoliubov transformation in 
which new fermion operators 7^ and 7^ are introduced 
as linear combination of Ck and ct 



i=i 



(34) 



7fc = u k c k - iv k 4_ k , 
7fe = u k4 + iv k c-k, 



(39) 
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where u k and v k are real numbers parametrized by 
Uk = sin ^ and V/. = cos % . This parametrization nat- 
urally arises as a consequence of the fermionic algebra 
{lk,ll>} = S kk >, {%, = {ik,lk<} = 0, with 5 kk , 
standing for the Kronecker delta symbol. Moreover, to 
recast the Hamiltonian in a diagonal form we define Ok 
by demanding that tan#ft = sin ka/(g — cos ka). There- 
fore, by expressing Hj in terms of Bogoliubov fcrmions 
and by imposing the trace invariance of the Hamiltonian, 
Eq. (J35J) becomes 



#7 = $> 



t 1 



(40) 



with Ek = 2Jyl + g 2 — 2g cos ka. Hamiltonian ([40]) is 
diagonal, with ground state given by the 7-fermion vac- 
uum. The procedure above also applies for the evalua- 
tion of the matrix elements of the reduced density op- 
erator given by Eq. (|12p , which amounts for the com- 
putation of the magnetization density G z and the two- 
point functions G a p. This can be achieved by using that 
G zz = G 2 z — G xx G yy [13] and by expressing the remaining 
correlation functions as 



M 



Gxx 

G z 
where 



/2vrg\ 2 . f2Trq\ 

COS I —jT- J V q + Sm I J U 9 V 9 



E 

q=-M 
M 

E 

= -7! 
I M 

T E (!-^ 2 )' 



q=-M 
M 



(2nq\ 2 . (2nq\ 
cos I — 1 u g - sm I — 1 



(41) 



l + . 9 2-2 3 cos(^) 



^l + ff 2-2 ff cos(^) 



(42) 



Hence, we exactly determine the two-spin reduced den- 
sity matrix. Classical and quantum correlations can then 
be directly obtained from Eqs. ([9|) and (fT0|) . By numer- 
ically computing the classical correlation in Eq. (|25p for 
nearest- neighbor spin pairs at sites i and i + 1, we can 
show that the maximization is achieved for any g by the 
choice 9 = ir/2 and cf> = 0. Then, the measurement that 
maximizes J(p : {-Bfc}) is given by {|+)(+|, |— )(— |}, with 
|+) and |— ) denoting the up and down spins in the x di- 
rection, namely, |±) = (| \) ± | X))j\[2. This numerical 
observation implies that w\ = 1, w% = W3 = 0. There- 
fore, Eq. (|25|) is ruled by the spin functions Ci = G^.'* +1 
and C4 = C5 = G\, i.e. 



where _fff,i„ is the binary entropy 

Hbinip) = -plogp - (1 - p) log (1 - p) 

and 

1 

Pi = 
P2 = 



(44) 




(45) 



We plot C(p) and Q(p) in Fig. [3] for a chain with 1024 
sites. Note that, for g = the system is a classical Ising 
chain, whose ground state is a doublet given by the vec- 
tors ! + +•••+) and I •••—). Therefore, the system 

is in the mixed state 



p = ^\ + + ---+)(+ + --- + \ + ^\ X 1, 

(46) 

with C(p) = I(p) = 1 and Q(p) — 0. On the other hand, 
in the limit g — > 00 the system is a paramagnet (vanish- 
ing magnetization in the x direction), with all spins in 
state I \). Therefore the system will be described by the 
density operator 



P = |TT • - ■ TXTT — T U 



(47) 



which is a pure separable state, containing neither clas- 
sical nor quantum correlations. 



Classical Correlation 
Quantum Correlation 




C(p) = Hun (Pi) - Hbin (P2) 



(43) 



FIG. 3: (Color online) Classical and quantum correlations for 
nearest-neighbor spins in the transverse field Ising model for 
a chain with 1024 sites. 

The QPT from ferromagnetic to paramagnetic state is 
a second-order QPT and occurs at g = 1. Signatures 
of this QPT can be found out by looking at the deriva- 
tives of either classical or quantum correlations. Indeed, 
the QPT can be identified as a pronounced minimum of 
the first derivative of the classical correlation, which is 
exhibited in Fig. [4] Note that the minimum logarithmi- 
cally diverges at g = 1 as the thermodynamic limit is 
approached (see inset of Fig. QI • In the case of quantum 
correlations, its first derivative shows an inflexion point 



0.5 1 1.5 2 

g 

FIG. 4: (Color online) First derivative of the classical correla- 
tion for nearest-neighbor spins with respect to g in the trans- 
verse field Ising chain for different lattice sizes L. The deriva- 
tive of C has a pronounced minimum at g m in, which tends 
to the QCP g = I as L — > oo. Inset: dC/dg taken at g m in 
exhibits a logarithmic divergence fitted by (dC/dg)\ = 
—0.29161 — 0.22471 log L. 

around g = 1, as displayed in Fig. [5] Indeed, by looking 
at its second derivative in Fig. the QPT is identified 
by a pronounced maximum, which shows quadratic loga- 
rithmic divergence at g = 1 as the thermodynamic limit 
is approached (see inset of Fig. [5]). 
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FIG. 5: (Color online) First derivative of the quantum cor- 
relation for nearest-neighbor spins with respect to g in the 
transverse field Ising chain for different lattice sizes L. In- 
set: dQ/dg presents an inflexion point that tends to the QCP 
g = 1 as L — > oo. 
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FIG. 6: (Color online) Second derivative of the quantum cor- 
relations for nearest-neighbor spins with respect to g in the 
transverse field Ising chain for different lattice sizes L. Ob- 
serve that the second derivative of Q has a pronounced max- 
imum at g m ax, which tends to the QCP g = 1 as L — > oo. 
Inset: d 2 C/dg 2 taken at g ma x exhibits a quadratic logarithmic 
divergence fitted by {d 2 Q/dg 2 )\ = 1.268-0.94712 logL+ 
0.15176 log 2 L. 



esting to observe that the spin functions Ci = G^ +1 and 
C4 = C5 = G l zl which govern the correlations in the Ising 
chain [see Eqs. (|43[) - (|4"5|) ]. exhibit a crossing at the QCP. 
This is shown in Fig. [7] for a chain with 1024 sites. 




FIG. 7: (Color online) Spin functions ci = G^ +1 and ca = G\ 
as a function of g. These functions govern the behavior of 
correlations, exhibiting a crossing exactly at the QCP. The 
plot is for a chain with L — 1024 sites. 



The behavior of the quantum discord is therefore 
rather different from the entanglement behavior, whose 
first derivative is already divergent at the QCP. Remark- 
ably, the scaling of pairwise entanglement derivative in 
this case (see e.g. Refs. is much closer to the scaling 

of the classical correlation derivative (as given by Fig. 
than that of the quantum correlation derivative (as given 
by Fig. [5]). As in the case of the XXZ model, it is inter- 



VI. THE LMG MODEL 

The discussion of correlations above can also be ap- 
plied in collective systems. As an illustration, we will 
consider here the LMG model 38], which describes a 
two-level Fermi system {|+),|— )}, with each level hav- 
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ing degeneracy The Hamiltonian for LMG model is 
given by 



H = X 



2 ( C +m C +™ C -m C -m) 

m— 1 
m,n— 1 



The operators ct TO and c_ m create a particle in the up- 
per and lower levels, respectively This Hamiltonian can 
be taken as describing an effective model for many-body 
systems, with one level just below the Fermi level and 
and the other level just above, with the level below be- 
ing filled with f2 particles [3^|. Alternatively, the LMG 
model can be seen as a one-dimensional ring of spin-1/2 
particles with infinite range interaction between pairs. 
Indeed, the Hamiltonian can be rewritten as 



(48) 



H - XS Z - — (Si 



(49) 



where S z = 



N 



and S x +iSy = 



12 ( c +m c +m — C_ m C. 

53rn=i c \m c -m ,39]. The system undergoes a second- 
order QPT at A = 1. As 57 — > oo, the ground state, as 
given by the Hartree-Fock (HF) approach, reads 



\HF) = Yl al 



(50) 



where we have introduced new levels labelled by and 1 
governed by the operators 



cos a c_ m + sin a c\ m , 



= — sin a c_ 



cos a c 



+ 771 ' 



(51) 



In Eq. (|51[) . a is a variational parameter to be adjusted 
in order to minimize energy, which is achieved according 
to the choice 



A < 1 
A > 1 



cos 2a = A, 
a = 0. 



(52) 



Despite being an approximation, the HF ground state 
provides the exact description of the critical point (for 
recent discussions of the exact spectrum of the LMG 
model, see Refs. [53, EU). The pairwise density oper- 
ator for general modes i = (+m) and j = (— n) is given 
by 



/ (MiMj) 



V o 



(MiNj 





(ch) 

(NtMj) 


t. 









(53) 



where Mk = 1 — and Nk = c[ck, with k — By 
evaluating the matrix elements of p for the HF ground 



state, we obtain 

(M +m M- n ) 
(M +m N_ n ) 
(N +m M_ n ) 
(N +m N_ n ) 

( c +m c -n) 
( c -ri c +m) 



sin 2 a cos 2 a (1 — 8 mn ) , 
cos 2 aS mn + cos 4 a (1 - 5 mn ) 
sin 2 a5 mn + sin 4 a (1 - <5 TOn ) 



= sin 2 a cos 2 a (1 
= sin a cos aS mn 
= sin a cos ad mn . 







(54) 



Note that Eq. (|55^) displays Z2 symmetry and, therefore, 
classical and quantum correlations can be computed by 
using Eq. ((23|) . Note also that, for m 7^ n, the density 
matrix is diagonal and the state is completely pairwise 
uncorrelated. On the other hand, for m — n, there is 
an equal amount of classical and quantum correlations 
between the modes. These correlations vanish for A > 1, 
which is the fully polarized state. The result is plotted 
in Fig. [5] We can then observe that the derivatives of 
both classical correlation and quantum discord exhibit a 
signature of the QPT (see inset of Fig. [H]) . These signa- 
tures are in agreement with the caracterizations in terms 
of entanglement [12, [H| and Fisher information [44[ . 




FIG. 8: (Color online) Classical and quantum correlations 
between modes (+m) and (—m) in the HF ground state of the 
Lipkin model. Inset: Derivative of correlations is nonanalytic 
at A = 1. 



VII. CONCLUSION 

In conclusion, we have investigated the behavior of 
pairwise correlations in general Z2 symmetric systems, 
splitting up their classical and quantum contributions. 
This allowed for the treatment of spin systems such as the 
XXZ model and the transverse field Ising chain, where 
the maximization required for the evaluation of the cor- 
relations has been analytically worked out. Moreover, we 
have identified the spin functions that govern the corre- 
lations and also discussed their behavior for finite size 
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chains. As a further application, we have used our ap- 
proach to investigate the case of a collective model, given 
by the LMG Hamiltonian. 

As it was shown, both classical correlation and quan- 
tum discord display signatures of the critical behavior of 
the system for the cases of first-order, second-order, and 
infinite-order QPTs. For first-order QPTs, as illustrated 
by the ferromagnetic point of the XXZ spin chain, we 
found that both classical correlation and quantum dis- 
cord display a jump at the critical point, which closely re- 
sembles the behavior of entanglement. For second-order 
QPTs, although both kind of correlations exhibit signa- 
tures of the QPTs, the derivatives of the quantum discord 
show a scaling that is rather different from the entangle- 
ment behavior in the transverse field Ising model (see, 
e.g., Refs. [HHI). It is remarkable that, for this model, 
the first derivative of pairwise entanglement with respect 
to the parameter that drives the QPT exhibits a scaling 
that is much closer to the scaling of the first derivative 
of the classical correlation. For infinite-order QPTs, as 
given by the antiferromagnetic point of the XXZ spin 



chain, classical correlation is a minimum at the QCP 
while quantum discord is a maximum. Whether or not 
this might be a general feature of infinite-order QPTs is 
still under analysis. 

Further investigations including correlations between 
blocks of particles and the effect of temperature may 
be interesting to establish a precise comparison between 
(classical and quantum) correlations and entanglement 
at QPTs. Moreover, dynamics in open quantum sys- 
tems p5l l46l |47| may also provide an interesting scenario 
for the discussion of the properties of the correlations and 
its implications for phase transitions. Such topics are left 
for a future research. 
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